shaPRS: Leveraging shared genetic effects across traits or ancestries improves accuracy of polygenic scores.

We present shaPRS, a method that leverages widespread pleiotropy between traits or shared genetic effects across ancestries, to improve the accuracy of polygenic scores. The method uses genome-wide summary statistics from two diseases or ancestries to improve the genetic effect estimate and standard error at SNPs where there is homogeneity of effect between the two datasets. When there is significant evidence of heterogeneity, the genetic effect from the disease or population closest to the target population is maintained. We show via simulation and a series of real-world examples that shaPRS substantially enhances the accuracy of polygenic risk scores (PRSs) for complex diseases and greatly improves PRS performance across ancestries. shaPRS is a PRS pre-processing method that is agnostic to the actual PRS generation method, and as a result, it can be integrated into existing PRS generation pipelines and continue to be applied as more performant PRS methods are developed over time.

Large-scale meta-genome-wide association study reveals common genetic factors linked to radiation-induced acute toxicities across cancer types.

BACKGROUND: This study was designed to identify common genetic susceptibility and shared genetic variants associated with acute radiation-induced toxicity across 4 cancer types (prostate, head and neck, breast, and lung). METHODS: A genome-wide association study meta-analysis was performed using 19 cohorts totaling 12 042 patients. Acute standardized total average toxicity (STATacute) was modelled using a generalized linear regression model for additive effect of genetic variants, adjusted for demographic and clinical covariates (rSTATacute). Linkage disequilibrium score regression estimated shared single-nucleotide variation (SNV-formerly SNP)-based heritability of rSTATacute in all patients and for each cancer type. RESULTS: Shared SNV-based heritability of STATacute among all cancer types was estimated at 10% (SE = 0.02) and was higher for prostate (17%, SE = 0.07), head and neck (27%, SE = 0.09), and breast (16%, SE = 0.09) cancers. We identified 130 suggestive associated SNVs with rSTATacute (5.0 × 10‒8 < P < 1.0 × 10‒5) across 25 genomic regions. rs142667902 showed the strongest association (effect allele A; effect size ‒0.17; P = 1.7 × 10‒7), which is located near DPPA4, encoding a protein involved in pluripotency in stem cells, which are essential for repair of radiation-induced tissue injury. Gene-set enrichment analysis identified 'RNA splicing via endonucleolytic cleavage and ligation' (P = 5.1 × 10‒6, P = .079 corrected) as the top gene set associated with rSTATacute among all patients. In silico gene expression analysis showed that the genes associated with rSTATacute were statistically significantly up-regulated in skin (not sun exposed P = .004 corrected; sun exposed P = .026 corrected). CONCLUSIONS: There is shared SNV-based heritability for acute radiation-induced toxicity across and within individual cancer sites. Future meta-genome-wide association studies among large radiation therapy patient cohorts are worthwhile to identify the common causal variants for acute radiotoxicity across cancer types.

Genome-wide association study of treatment-related toxicity two years following radiotherapy for breast cancer.

BACKGROUND AND PURPOSE: Up to a quarter of breast cancer patients treated by surgery and radiotherapy experience clinically significant toxicity. If patients at high risk of adverse effects could be identified at diagnosis, their treatment could be tailored accordingly. This study was designed to identify common single nucleotide polymorphisms (SNPs) associated with toxicity two years following whole breast radiotherapy. MATERIALS AND METHODS: A genome-wide association study (GWAS) was performed in 1,640 breast cancer patients with complete SNP, clinical, treatment and toxicity data, recruited across 18 European and US centres into the prospective REQUITE cohort study. Toxicity data (CTCAE v4.0) were collected at baseline, end of radiotherapy, and annual follow-up. A total of 7,097,340 SNPs were tested for association with the residuals of toxicity endpoints, adjusted for clinical, treatment co-variates and population substructure. RESULTS: Quantile-quantile plots showed more associations with toxicity above the p < 5 × 10-5 level than expected by chance. Eight SNPs reached genome-wide significance. Nipple retraction grade ≥ 2 was associated with the rs188287402 variant (p = 2.80 × 10-8), breast oedema grade ≥ 2 with rs12657177 (p = 1.12 × 10-10), rs75912034 (p = 1.12 × 10-10), rs145328458 (p = 1.06 × 10-9) and rs61966612 (p = 1.23 × 10-9), induration grade ≥ 2 with rs77311050 (p = 2.54 × 10-8) and rs34063419 (p = 1.21 × 10-8), and arm lymphoedema grade ≥ 1 with rs643644 (p = 3.54 × 10-8). Heritability estimates across significant endpoints ranged from 25% to 39%. Our study did not replicate previously reported SNPs associated with breast radiation toxicity at the pre-specified significance level. CONCLUSIONS: This GWAS for long-term breast radiation toxicity provides further evidence for significant association of common SNPs with distinct toxicity endpoints.

Treatment time and circadian genotype interact to influence radiotherapy side-effects. A prospective European validation study using the REQUITE cohort.

BACKGROUND: Circadian rhythm impacts broad biological processes, including response to cancer treatment. Evidence conflicts on whether treatment time affects risk of radiotherapy side-effects, likely because of differing time analyses and target tissues. We previously showed interactive effects of time and genotypes of circadian genes on late toxicity after breast radiotherapy and aimed to validate those results in a multi-centre cohort. METHODS: Clinical and genotype data from 1690 REQUITE breast cancer patients were used with erythema (acute; n=340) and breast atrophy (two years post-radiotherapy; n=514) as primary endpoints. Local datetimes per fraction were converted into solar times as predictors. Genetic chronotype markers were included in logistic regressions to identify primary endpoint predictors. FINDINGS: Significant predictors for erythema included BMI, radiation dose and PER3 genotype (OR 1.27(95%CI 1.03-1.56); P < 0.03). Effect of treatment time effect on acute toxicity was inconclusive, with no interaction between time and genotype. For late toxicity (breast atrophy), predictors included BMI, radiation dose, surgery type, treatment time and SNPs in CLOCK (OR 0.62 (95%CI 0.4-0.9); P < 0.01), PER3 (OR 0.65 (95%CI 0.44-0.97); P < 0.04) and RASD1 (OR 0.56 (95%CI 0.35-0.89); P < 0.02). There was a statistically significant interaction between time and genotypes of circadian rhythm genes (CLOCK OR 1.13 (95%CI 1.03-1.23), P < 0.01; PER3 OR 1.1 (95%CI 1.01-1.2), P < 0.04; RASD1 OR 1.15 (95%CI 1.04-1.28), P < 0.008), with peak time for toxicity determined by genotype. INTERPRETATION: Late atrophy can be mitigated by selecting optimal treatment time according to circadian genotypes (e.g. treat PER3 rs2087947C/C genotypes in mornings; T/T in afternoons). We predict triple-homozygous patients (14%) reduce chance of atrophy from 70% to 33% by treating in mornings as opposed to mid-afternoon. Future clinical trials could stratify patients treated at optimal times compared to those scheduled normally. FUNDING: EU-FP7.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

No Association Between Polygenic Risk Scores for Cancer and Development of Radiation Therapy Toxicity.

PURPOSE: Our aim was to test whether updated polygenic risk scores (PRS) for susceptibility to cancer affect risk of radiation therapy toxicity. METHODS AND MATERIALS: Analyses included 9,717 patients with breast (n=3,078), prostate (n=5,748) or lung (n=891) cancer from Radiogenomics and REQUITE Consortia cohorts. Patients underwent potentially curative radiation therapy and were assessed prospectively for toxicity. Germline genotyping involved genome-wide single nucleotide polymorphism (SNP) arrays with nontyped SNPs imputed. PRS for each cancer were generated by summing literature-identified cancer susceptibility risk alleles: 352 breast, 136 prostate, and 24 lung. Weighted PRS were generated using log odds ratio (ORs) for cancer susceptibility. Standardized total average toxicity (STAT) scores at 2 and 5 years (breast, prostate) or 6 to 12 months (lung) quantified toxicity. Primary analysis tested late STAT, secondary analyses investigated acute STAT, and individual endpoints and SNPs using multivariable regression. RESULTS: Increasing PRS did not increase risk of late toxicity in patients with breast (OR, 1.000; 95% confidence interval [CI], 0.997-1.002), prostate (OR, 0.99; 95% CI, 0.98-1.00; weighted PRS OR, 0.93; 95% CI, 0.83-1.03), or lung (OR, 0.93; 95% CI, 0.87-1.00; weighted PRS OR, 0.68; 95% CI, 0.45-1.03) cancer. Similar results were seen for acute toxicity. Secondary analyses identified rs138944387 associated with breast pain (OR, 3.05; 95% CI, 1.86-5.01; P = 1.09 × 10-5) and rs17513613 with breast edema (OR, 0.94; 95% CI, 0.92-0.97; P = 1.08 × 10-5). CONCLUSIONS: Patients with increased polygenic predisposition to breast, prostate, or lung cancer can safely undergo radiation therapy with no anticipated excess toxicity risk. Some individual SNPs increase the likelihood of a specific toxicity endpoint, warranting validation in independent cohorts and functional studies to elucidate biologic mechanisms.

Genome-wide association study and functional follow-up identify 14q12 as a candidate risk locus for cervical cancer.

Cervical cancer is among the leading causes of cancer-related death in females worldwide. Infection by human papillomavirus (HPV) is an established risk factor for cancer development. However, genetic factors contributing to disease risk remain largely unknown. We report on a genome-wide association study (GWAS) on 375 German cervical cancer patients and 866 healthy controls, followed by a replication study comprising 658 patients with invasive cervical cancer, 1361 with cervical dysplasia and 841 healthy controls. Functional validation was performed for the top GWAS variant on chromosome 14q12 (rs225902, close to PRKD1). After bioinformatic annotation and in silico predictions, we performed transcript analysis in a cervical tissue series of 317 samples and demonstrate rs225902 as an expression quantitative trait locus (eQTL) for FOXG1 and two tightly co-regulated long non-coding RNAs at this genomic region, CTD-2251F13 (lnc-PRKD1-1) and CTD-2503I6 (lnc-FOXG1-6). We also show allele-specific effects of the 14q12 variants via luciferase assays. We propose a combined effect of genotype, HPV status and gene expression at this locus on cervical cancer progression. Taken together, this work uncovers a potential candidate locus with regulatory functions and contributes to the understanding of genetic susceptibility to cervical cancer.

A genome-wide association study of radiotherapy induced toxicity in head and neck cancer patients identifies a susceptibility locus associated with mucositis.

PURPOSE: A two-stage genome-wide association study was carried out in head and neck cancer (HNC) patients aiming to identify genetic variants associated with either specific radiotherapy-induced (RT) toxicity endpoints or a general proneness to develop toxicity after RT. MATERIALS AND METHODS: The analysis included 1780 HNC patients treated with primary RT for laryngeal or oro/hypopharyngeal cancers. In a non-hypothesis-driven explorative discovery study, associations were tested in 1183 patients treated within The Danish Head and Neck Cancer Group. Significant associations were later tested in an independent Dutch cohort of 597 HNC patients and if replicated, summary data obtained from discovery and replication studies were meta-analysed. Further validation of significantly replicated findings was pursued in an Asian cohort of 235 HNC patients with nasopharynx as the primary tumour site. RESULTS: We found and replicated a significant association between a locus on chromosome 5 and mucositis with a pooled OR for rs1131769*C in meta-analysis = 1.95 (95% CI 1.48-2.41; ppooled = 4.34 × 10-16). CONCLUSION: This first exploratory GWAS in European cohorts of HNC patients identified and replicated a risk locus for mucositis. A larger Meta-GWAS to identify further risk variants for RT-induced toxicity in HNC patients is warranted.

Functional annotation of the 2q35 breast cancer risk locus implicates a structural variant in influencing activity of a long-range enhancer element.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

Development of a method for generating SNP interaction-aware polygenic risk scores for radiotherapy toxicity.

AIM: To identify the effect of single nucleotide polymorphism (SNP) interactions on the risk of toxicity following radiotherapy (RT) for prostate cancer (PCa) and propose a new method for polygenic risk score incorporating SNP-SNP interactions (PRSi). MATERIALS AND METHODS: Analysis included the REQUITE PCa cohort that received external beam RT and was followed for 2 years. Late toxicity endpoints were: rectal bleeding, urinary frequency, haematuria, nocturia, decreased urinary stream. Among 43 literature-identified SNPs, the 30% most strongly associated with each toxicity were tested. SNP-SNP combinations (named SNP-allele sets) seen in ≥10% of the cohort were condensed into risk (RS) and protection (PS) scores, respectively indicating increased or decreased toxicity risk. Performance of RS and PS was evaluated by logistic regression. RS and PS were then combined into a single PRSi evaluated by area under the receiver operating characteristic curve (AUC). RESULTS: Among 1,387 analysed patients, toxicity rates were 11.7% (rectal bleeding), 4.0% (urinary frequency), 5.5% (haematuria), 7.8% (nocturia) and 17.1% (decreased urinary stream). RS and PS combined 8 to 15 different SNP-allele sets, depending on the toxicity endpoint. Distributions of PRSi differed significantly in patients with/without toxicity with AUCs ranging from 0.61 to 0.78. PRSi was better than the classical summed PRS, particularly for the urinary frequency, haematuria and decreased urinary stream endpoints. CONCLUSIONS: Our method incorporates SNP-SNP interactions when calculating PRS for radiotherapy toxicity. Our approach is better than classical summation in discriminating patients with toxicity and should enable incorporating genetic information to improve normal tissue complication probability models.

Breast Cancer Risk Genes - Association Analysis in More than 113,000 Women.

BACKGROUND: Genetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking. METHODS: We used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity. RESULTS: Protein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants. CONCLUSIONS: The results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.).

A Deep Learning Approach Validates Genetic Risk Factors for Late Toxicity After Prostate Cancer Radiotherapy in a REQUITE Multi-National Cohort.

Background: REQUITE (validating pREdictive models and biomarkers of radiotherapy toxicity to reduce side effects and improve QUalITy of lifE in cancer survivors) is an international prospective cohort study. The purpose of this project was to analyse a cohort of patients recruited into REQUITE using a deep learning algorithm to identify patient-specific features associated with the development of toxicity, and test the approach by attempting to validate previously published genetic risk factors. Methods: The study involved REQUITE prostate cancer patients treated with external beam radiotherapy who had complete 2-year follow-up. We used five separate late toxicity endpoints: ≥grade 1 late rectal bleeding, ≥grade 2 urinary frequency, ≥grade 1 haematuria, ≥ grade 2 nocturia, ≥ grade 1 decreased urinary stream. Forty-three single nucleotide polymorphisms (SNPs) already reported in the literature to be associated with the toxicity endpoints were included in the analysis. No SNP had been studied before in the REQUITE cohort. Deep Sparse AutoEncoders (DSAE) were trained to recognize features (SNPs) identifying patients with no toxicity and tested on a different independent mixed population including patients without and with toxicity. Results: One thousand, four hundred and one patients were included, and toxicity rates were: rectal bleeding 11.7%, urinary frequency 4%, haematuria 5.5%, nocturia 7.8%, decreased urinary stream 17.1%. Twenty-four of the 43 SNPs that were associated with the toxicity endpoints were validated as identifying patients with toxicity. Twenty of the 24 SNPs were associated with the same toxicity endpoint as reported in the literature: 9 SNPs for urinary symptoms and 11 SNPs for overall toxicity. The other 4 SNPs were associated with a different endpoint. Conclusion: Deep learning algorithms can validate SNPs associated with toxicity after radiotherapy for prostate cancer. The method should be studied further to identify polygenic SNP risk signatures for radiotherapy toxicity. The signatures could then be included in integrated normal tissue complication probability models and tested for their ability to personalize radiotherapy treatment planning.

Somatic mosaicism and common genetic variation contribute to the risk of very-early-onset inflammatory bowel disease.

Very-early-onset inflammatory bowel disease (VEO-IBD) is a heterogeneous phenotype associated with a spectrum of rare Mendelian disorders. Here, we perform whole-exome-sequencing and genome-wide genotyping in 145 patients (median age-at-diagnosis of 3.5 years), in whom no Mendelian disorders were clinically suspected. In five patients we detect a primary immunodeficiency or enteropathy, with clinical consequences (XIAP, CYBA, SH2D1A, PCSK1). We also present a case study of a VEO-IBD patient with a mosaic de novo, pathogenic allele in CYBB. The mutation is present in ~70% of phagocytes and sufficient to result in defective bacterial handling but not life-threatening infections. Finally, we show that VEO-IBD patients have, on average, higher IBD polygenic risk scores than population controls (99 patients and 18,780 controls; P < 4 × 10-10), and replicate this finding in an independent cohort of VEO-IBD cases and controls (117 patients and 2,603 controls; P < 5 × 10-10). This discovery indicates that a polygenic component operates in VEO-IBD pathogenesis.

Novel CYP4F22 mutations associated with autosomal recessive congenital ichthyosis (ARCI). Study of the CYP4F22 c.1303C>T founder mutation.

Mutations in CYP4F22 cause autosomal recessive congenital ichthyosis (ARCI). However, less than 10% of all ARCI patients carry a mutation in CYP4F22. In order to identify the molecular basis of ARCI among our patients (a cohort of ninety-two Spanish individuals) we performed a mutational analysis using direct Sanger sequencing in combination with a multigene targeted NGS panel. From these, eight ARCI families (three of them with Moroccan origin) were found to carry five different CYP4F22 mutations, of which two were novel. Computational analysis showed that the mutations found were present in highly conserved residues of the protein and may affect its structure and function. Seven of the eight families were carriers of a highly recurrent CYP4F22 variant, c.1303C>T; p.(His435Tyr). A 12Mb haplotype was reconstructed in all c.1303C>T carriers by genotyping ten microsatellite markers flanking the CYP4F22 gene. A prevalent 2.52Mb haplotype was observed among Spanish carrier patients suggesting a recent common ancestor. A smaller core haplotype of 1.2Mb was shared by Spanish and Moroccan families. Different approaches were applied to estimate the time to the most recent common ancestor (TMRCA) of carrier patients with Spanish origin. The age of the mutation was calculated by using DMLE and BDMC2. The algorithms estimated that the c.1303C>T variant arose approximately 2925 to 4925 years ago, while Spanish carrier families derived from a common ancestor who lived in the XIII century. The present study reports five CYP4F22 mutations, two of them novel, increasing the number of CYP4F22 mutations currently listed. Additionally, our results suggest that the recurrent c.1303C>T change has a founder effect in Spanish population and c.1303C>T carrier families originated from a single ancestor with probable African ancestry.

Candidate Causal Variants at the 8p12 Breast Cancer Risk Locus Regulate DUSP4.

Genome-wide association studies have revealed a locus at 8p12 that is associated with breast cancer risk. Fine-mapping of this locus identified 16 candidate causal variants (CCVs). However, as these variants are intergenic, their function is unclear. To map chromatin looping from this risk locus to a previously identified candidate target gene, DUSP4, we performed chromatin conformation capture analyses in normal and tumoural breast cell lines. We identified putative regulatory elements, containing CCVs, which looped to the DUSP4 promoter region. Using reporter gene assays, we found that the risk allele of CCV rs7461885 reduced the activity of a DUSP4 enhancer element, consistent with the function of DUSP4 as a tumour suppressor gene. Furthermore, the risk allele of CCV rs12155535, located in another DUSP4 enhancer element, was negatively correlated with looping of this element to the DUSP4 promoter region, suggesting that this allele would be associated with reduced expression. These findings provide the first evidence that CCV risk alleles downregulate DUSP4 expression, suggesting that this gene is a regulatory target of the 8p12 breast cancer risk locus.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes.

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Chromatin interactome mapping at 139 independent breast cancer risk signals.

BACKGROUND: Genome-wide association studies have identified 196 high confidence independent signals associated with breast cancer susceptibility. Variants within these signals frequently fall in distal regulatory DNA elements that control gene expression. RESULTS: We designed a Capture Hi-C array to enrich for chromatin interactions between the credible causal variants and target genes in six human mammary epithelial and breast cancer cell lines. We show that interacting regions are enriched for open chromatin, histone marks for active enhancers, and transcription factors relevant to breast biology. We exploit this comprehensive resource to identify candidate target genes at 139 independent breast cancer risk signals and explore the functional mechanism underlying altered risk at the 12q24 risk region. CONCLUSIONS: Our results demonstrate the power of combining genetics, computational genomics, and molecular studies to rationalize the identification of key variants and candidate target genes at breast cancer GWAS signals.

Non-coding RNAs underlie genetic predisposition to breast cancer.

BACKGROUND: Genetic variants identified through genome-wide association studies (GWAS) are predominantly non-coding and typically attributed to altered regulatory elements such as enhancers and promoters. However, the contribution of non-coding RNAs to complex traits is not clear. RESULTS: Using targeted RNA sequencing, we systematically annotated multi-exonic non-coding RNA (mencRNA) genes transcribed from 1.5-Mb intervals surrounding 139 breast cancer GWAS signals and assessed their contribution to breast cancer risk. We identify more than 4000 mencRNA genes and show their expression distinguishes normal breast tissue from tumors and different breast cancer subtypes. Importantly, breast cancer risk variants, identified through genetic fine-mapping, are significantly enriched in mencRNA exons, but not the promoters or introns. eQTL analyses identify mencRNAs whose expression is associated with risk variants. Furthermore, chromatin interaction data identify hundreds of mencRNA promoters that loop to regions that contain breast cancer risk variants. CONCLUSIONS: We have compiled the largest catalog of breast cancer-associated mencRNAs to date and provide evidence that modulation of mencRNAs by GWAS variants may provide an alternative mechanism underlying complex traits.

Radiogenomics Consortium Genome-Wide Association Study Meta-Analysis of Late Toxicity After Prostate Cancer Radiotherapy.

BACKGROUND: A total of 10%-20% of patients develop long-term toxicity following radiotherapy for prostate cancer. Identification of common genetic variants associated with susceptibility to radiotoxicity might improve risk prediction and inform functional mechanistic studies. METHODS: We conducted an individual patient data meta-analysis of six genome-wide association studies (n = 3871) in men of European ancestry who underwent radiotherapy for prostate cancer. Radiotoxicities (increased urinary frequency, decreased urinary stream, hematuria, rectal bleeding) were graded prospectively. We used grouped relative risk models to test associations with approximately 6 million genotyped or imputed variants (time to first grade 2 or higher toxicity event). Variants with two-sided Pmeta less than 5 × 10-8 were considered statistically significant. Bayesian false discovery probability provided an additional measure of confidence. Statistically significant variants were evaluated in three Japanese cohorts (n = 962). All statistical tests were two-sided. RESULTS: Meta-analysis of the European ancestry cohorts identified three genomic signals: single nucleotide polymorphism rs17055178 with rectal bleeding (Pmeta = 6.2 × 10-10), rs10969913 with decreased urinary stream (Pmeta = 2.9 × 10-10), and rs11122573 with hematuria (Pmeta = 1.8 × 10-8). Fine-scale mapping of these three regions was used to identify another independent signal (rs147121532) associated with hematuria (Pconditional = 4.7 × 10-6). Credible causal variants at these four signals lie in gene-regulatory regions, some modulating expression of nearby genes. Previously identified variants showed consistent associations (rs17599026 with increased urinary frequency, rs7720298 with decreased urinary stream, rs1801516 with overall toxicity) in new cohorts. rs10969913 and rs17599026 had similar effects in the photon-treated Japanese cohorts. CONCLUSIONS: This study increases the understanding of the architecture of common genetic variants affecting radiotoxicity, points to novel radio-pathogenic mechanisms, and develops risk models for testing in clinical studies. Further multinational radiogenomics studies in larger cohorts are worthwhile.